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We determine the strong coupling constant a s and its energy dependence from the pr dependence 
of the inclusive jet cross section in pp collisions at \fs — 1.96 TeV. The strong coupling constant 
is determined over the transverse momentum range 50 < pr < 145 GeV. Using perturbative QCD 
calculations to order 0(« 3 ) combined with 0(a 4 ) contributions from threshold corrections, we obtain 
a s (Mz) = 0.11611q Q048- This is the most precise result obtained at a hadron-hadron collider. 

PACS numbers: 13.87.-a, 12.38.Qk, 13.87.Ce 



Asymptotic freedom, the fact that the strong force be- 
tween quarks and gluons keeps getting weaker when it 
is probed at increasingly small distances, is a remark- 
able property of quantum chromodynamics (QCD). This 
property is reflected by the renormalization group equa- 
tion (RGE) prediction for the dependence of the strong 
coupling constant a s on the renormalization scale \x r 
and therefore on the momentum transfer. Experimental 
tests of asymptotic freedom require precise determina- 
tions of a s (ii r ) over a large range of momentum trans- 



fer. Frequently, a s has been determined using produc- 
tion rates of hadronic jets in either e + e~ annihilation 
or in deep- inelastic ep scattering (DIS) @. So far there 
exists only a single a s result from inclusive jet produc- 
tion in hadron-hadron collisions. The CDF Collabora- 
tion determined a s from the inclusive jet cross section 
in pp collisions at y/s — 1.8 TeV obtaining a s {M_z) = 
0.1178i8;8g||(exp.)±o:ooS(scalfi) ± 0.0059(PDF) 0. 

In this article we determine a s and its dependence 
on the momentum transfer using the published measure- 
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ment of the inclusive jet cross section [ll|, [T3| with the DO 
detector [l3| at the Fermilab Tevatron Collider in pp col- 
lisions at y/s = 1.96 TeV. The inclusive jet cross section 
d 2 a j e t /dpTd\y\ was measured using the Run II iterative 
midpoint cone algorithm (l4j with a cone radius of 0.7 in 
rapidity, y, and azimuthal angle. Rapidity is related to 
the polar scattering angle 9 with respect to the beam axis 
by y = 0.5 In [(1 + /? cos 0)/(l - /3 cos 9)] with /3 = 
The measurement comprises 110 data points corrected to 
the particle level [H[ and presented as a function of the 
momentum component transverse to the beam direction, 
Pt, for pt > 50 GeV in six regions of \y\ for < \y\ < 2.4. 

The ingredients of perturbative QCD (pQCD) calcu- 
lations in hadron collisions are a s , the perturbative co- 
efficients c n (in the n-th power of a s ), and the parton 
distribution functions (PDFs). Conceptually, PDFs de- 
pend only on the hadron momentum fraction x carried 
by the parton and on the factorization scale [if. In prac- 
tice, PDFs are determined from measurements of ob- 
servables which depend on a s . Therefore resulting PDF 
parametrizations depend on the assumption for a s made 
in the extraction procedure. For all precise phenomenol- 
ogy, this implicit a s dependence must be taken into ac- 
count consistently. The pQCD prediction for the inclu- 
sive jet cross section can therefore be written as 

CTport(a S ) = (^2 a7 s C nj ® fl( a s) ® M^s) > (1) 

where the sum runs over all powers n of a s which con- 
tribute to the calculation (n = 2,3,4 in this analysis, see 
below). The fi^ are the PDFs of the initial state hadrons 
and the "(g)" sign denotes the convolution over the mo- 
mentum fractions x±, xi of the hadrons. Since the RGE 
uniquely relates the value of a a (/j, r ) at any scale fx r to the 
value of a s (M z ), all equations can be expressed in terms 
of a s (Mz). The total theory prediction for inclusive jet 
production is given by the pQCD result in (JTJ) multiplied 
by a correction factor for nonperturbative effects 

— ^pcrt 

(a s (M z )) ■ 

^nonpert ■ (2) 

The factor c nonpcr t includes corrections due to hadroniza- 
tion and the underlying event which have been estimated 
in Ref. PH using PYTHIA pj] with CTEQ6.5 PDFs [13], 
tune QW [Hj], and a s (M z ) = 0.118. The hadroniza- 
tion (underlying event) corrections vary between -15% 
(+30%) to -3% (+6%), for p T = 50 to 600 GeV Q3]. 

The perturbative results are the sum of a full calcula- 
tion to 0(al) [next-to-leading order (NLO)], combined 
with the 0(a A s ) (2-loop) terms from threshold correc- 
tions . Adding the 2-loop threshold corrections leads 
to a significant reduction in the fi r and \x f dependence of 
the calculation. The theory calculations are performed 
in the MS scheme [2(| for five active quark flavors using 
the next-to-next-to-leading logarithmic (3-loop) approx- 
imation of the RGE 0, The PDFs are taken from 
the MSTW2008 next-to-next-to-leading order (NNLO) 
parametrizations [23l . [24| and fi r and \i / are both chosen 



equal to the jet pr- The calculations use fastnlo [25| 
based on NLOJET++ [26j, [27J and on code from the au- 
thors of Ref. [H. 

In this analysis, the value of a s is determined from 
sets of inclusive jet cross section data points by minimiz- 
ing the x 2 function between data and the theory result 
using MINUIT [28[. Where appropriate, the a s (M z ) 
result will be evolved to the scale px using the 3-loop 
solution of the RGE, providing a result for a s (pT)- All 
correlated experimental and theoretical uncertainties are 
treated in the Hessian approach |29j], except for the /i r j 
dependence (see below). The central a s (M z ) result is 
obtained by minimizing x 2 with respect to a s (M z ) and 
the nuisance parameters for the correlated uncertainties. 
By scanning x 2 as a function of a s (M z ), the uncertain- 
ties are obtained from the a s (M z ) values for which x 2 is 
increased by 1 with respect to the minimum value. 

To determine a s according to this procedure, knowl- 
edge of a per t(cts(M z )) is required as a continuous func- 
tion of a s (M z ), over a a s (M z ) range which covers the 
possible fit results and their uncertainties. This can be 
achieved based on a series of PDFs obtained under the 
same conditions but for different values of a s (M z ) us- 
ing interpolation in a s (M z ). Some recent PDF analy- 
ses have applied this strategy and their results are docu- 
mented for different values of a s {M z ). The MSTW2008 
NNLO (NLO) PDF parametrizations [H, [HI are pre- 
sented for 21 a s (M z ) values in the range 0.107 — 0.127 
(0.110 - 0.130) in steps of 0.001 and the CTEQ6.6 re- 
sults [3(j are available for five values of a s (M z ) = 
0.112,0.114,0.118,0.122,0.125. Because of the wide 
range in a s (M z ) covered by the MSTW2008 PDFs and 
the fine and equidistant spacing in a s (M z ), we use cu- 
bic spline interpolation to obtain a smooth parametriza- 
tion for the a s (M z ) dependence of the cross section for 
0.108 < a s (M z ) < 0.126 (0.111 < a s (M z ) < 0.129) 
for the NNLO (NLO) PDFs. This range is sufficient 
to cover our central values and the uncertainties. The 
MSTW2008 analysis includes data sets that have not 
yet been included in other global PDF analyses (DIS 
jet data from HERA and recent CCFR/NuTeV dimuon 
data); the results are available in NNLO accuracy which 
is adequate when including the 0(a A s ) contributions from 
threshold corrections in the cross section calculation. 
The CTEQ6.6 PDF parametrizations are available up 
to NLO, for five a s (M z ) values, and for a more limited 
range in a s (M z ) as compared to MSTW2008. Therefore 
the MSTW2008 PDFs are used to obtain the main re- 
sults for this analysis while the CTEQ6.6 PDFs are used 
for comparison. 

Care must be taken in phenomenological analyses if 
the observable under study was already used to provide 
significant constraints on the PDFs as this introduces 
correlations of experimental and PDF uncertainties, and 
it may affect the sensitivity to possible new physics sig- 
nals. Both aspects are relevant in this a s determination 
since the DO inclusive jet data under study is included in 
the MSTW2008 PDF analysis. Since the correlation of 
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experimental and PDF uncertainties is not documented, 
it can not be taken into account when using the PDFs 
to extract a s (Mz) from the jet data. As a consequence, 
we must avoid using those jet cross section data points 
which have provided strong PDF constraints. While the 
quark PDFs are constrained by precision structure func- 
tion data, the only direct source of information on the 
high x gluon PDF comes currently from Tevatron inclu- 
sive jet data. The impact of Tevatron jet data on the 
gluon density is documented in Ref. [HJ in Figs. 51-53. 
Figure 51 shows that excluding the Tevatron jet data 
starts to affect the gluon density at x > 0.2 — 0.3, while 
for x < 0.25 the difference in the gluon density with and 
without Tevatron jet data is less than 5%. Figure 53 
shows that x < 0.3 is the region in which the gluon re- 
sults for MSTW2008 and CTEQ6.6 are very close. We 
conclude that for momentum fractions x < 0.2 — 0.3 the 
Tevatron jet data do not have a significant impact on the 
gluon density, and therefore we can neglect correlations 
between PDF and experimental uncertainties for these 
data. Based on this constraint we select below those in- 
clusive jet data points from which we extract a s . 

The Tevatron jet data (which access px above 
500 GeV) are probing momentum transfers at which a s 
has not yet been probed in other experiments. There- 
fore we can not rule out deviations in the running of 
a s at large momentum due to possible new physics con- 
tributions to the RGE. Since such modifications of the 
RGE are not taken into account in the PDF determina- 
tions, these effects would effectively be absorbed into the 
PDFs. By construction, using such PDFs to extract a s 
could seemingly confirm the RGE expectations, even in 
the presence of new physics contributions to the RGE. 
For a consistent a s determination we would therefore ex- 
clude high pr data in the region where the RGE has 
not yet been successfully tested which is the region of 
Pt > 200 GeV 0. However, those data are already re- 
moved by the restriction to x < 0.2 — 0.3, so no additional 
requirement is needed to account for this. 

In 2 — > 2 processes, given the rapidities and px of the 
two jets, one can compute the momentum fractions X\ 
and Xi carried by the initial partons. The inclusive jet 
cross section at given pt and \y\ is, however, integrated 
over all additional jets in an event, so the rapidity of the 
other jet and therefore the full event kinematics, includ- 
ing X\ and X2 , are not known. The value of the larger mo- 
mentum fraction x max — max(a:i,X2) can be computed 
only under an assumption for the rapidity of the unob- 
served jet. For each inclusive jet (j>t, \y\) bin we define 



the variable x 



x T ■ (e 



\v\ 



l)/2 where xt = 2pt/\/s, 



Pt is taken at the bin center, and |y| at the lower bound- 
ary of the \y\ bin. This variable x corresponds to x ma x 
for the case that the unobserved jet was produced at 
y = 0. In the pQCD calculation, for a given inclusive 
jet (pt, \y\) bin the distribution of x max — max(:ri, £2) 
always has a peak plus a tail towards high a; maa; val- 
ues. Although the variable x does not represent the peak 
position of the x max distribution, it is correlated with 



0.2 



Q. 

^0 0.15 



0.1 

0.14 
0.12 
0.1 



N 



oc s (p T ) from inclusive jet cross section 
in hadron-induced processes 




B a s (M z )=0.116i:§X 
(D0 combined fit) 



4 I* 1 ^ } j {HHHH 



10 10 
Pt (GeV) 

FIG. 1: The results for a s (pT) (top) and a s (Mz) (bottom). 
The DO results are based on 22 selected data points which 
have been grouped to produce the 9 data points shown. For 
comparison, results from HERA DIS jet data have been in- 
cluded and also the RGE prediction for the combined DO fit 
result and its uncertainty (line and band) . All data points are 
shown with their total uncertainties. 



that distribution. The requirement x < 0.15 removes all 
data points for which more than half of the cross sec- 
tion is produced at x max > 0.25. This leaves 22 (out of 
110) data points for the a s analysis with pt < 145 GeV 
for < \y\ < 0.4, p T < 120 GeV for 0.4 < \y\ < 0.8, 
Pt < 90 GeV for 0.8 < \y\ < 1.2, and p T < 70 GeV for 
1.2 < \y\ < 1.6. Although this selection criterion is well- 
motivated, the specific choices of the variable x and the 
requirement x < 0.15 are somewhat arbitrary. We have 
therefore studied variations of the selection requirement 
in the range x < 0.10 — 0.17 and other choices for the defi- 
nition of x (for example assuming that the unobserved jet 
has 7/2 = =1= I y I ) , and, we find that the a s results are sta- 
ble within 1%. We conclude that the choice of x < 0.15 
restricts the jet data to those points which receive no sig- 
nificant contributions from x max > 0.25. For these data 
points, experimental and PDF uncertainties are treated 
as being uncorrelated. 

In the a s determination, we consider the uncorre- 
lated experimental uncertainties and all 23 sources of 
correlated experimental uncertainties as documented in 
Refs. [ll|, [T2] ■ The nonperturbative corrections are di- 
vided into hadronization and underlying event effects. 
The uncertainty for each is taken to be half the size of 
the corresponding effect. PDF uncertainties are com- 
puted using the 20 68% C.L. uncertainty eigenvectors 
as provided by MSTW2008 Q. The uncertainties in 
the pQCD calculation due to uncalculated higher order 
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TABLE I: Central values and uncertainties due to different sources for the nine a a {pr) results and for the combined a s (A4z) 
result (bottom). All uncertainties are multiplied by a factor of 10 3 . 



Pt range 


No. of 


Pt 




Total 


Experimental 


Experimental 


Nonperturb. 


PDF 




(GeV) 


data points 


(GeV) 


a s (pr) 


uncertainty 


uncorrelated 


correlated 


correction 


uncertainty 


variatk 


50-60 


4 


54.5 


0.1229 


+7.6 
— 7 7 


±0.4 


+4.8 
—4 9 


+ 5.8 

— 5 6 


+0.4 
— (> 


+ 1.0 

— 1 9 


60-70 


4 


64.5 


0.1204 


+6.2 

-6.3 


±0.3 


+4.1 
-4.3 


+4.5 
-4.3 


+0.6 
-0.5 


+ 1.3 
-1.5 


70-80 


3 


74.5 


0.1184 


+5.6 

-5.6 


±0.3 


+3.8 
-3.9 


+4.0 

-3.9 


+0.6 
-0.6 


+ 1.0 
-0.9 


80-90 


3 


84.5 


0.1163 


+5.1 
-5.1 


±0.3 


+3.6 
-3.7 


+3.5 

-3.5 


+0.7 
-0.7 


+0.9 
-0.6 


90-100 


2 


94.5 


0.1142 


+5.1 

-4.9 


±0.3 


+3.5 


+3.5 

-3.3 


+0.8 
-0.8 


+ 1.1 
-0.6 


100-110 


2 


104.5 


0.1131 


+4.7 
-4.7 


±0.2 


+1:1 

-3.5 


+3.1 
-3.0 


+0.8 
-0.8 


+ 1.1 
-0.6 


110-120 


2 


114.5 


0.1121 


+4.2 
-4.4 


±0.2 


+3.1 
-3.3 


+2.5 

-2.7 


+0.7 
-0.8 


+ 1.2 
-0.7 


120-130 


1 


124.5 


0.1102 


+4.4 
-4.4 


±0.2 


+3.2 
-3.4 


+2.6 

-2.6 


+0.9 
-0.9 


+ 1.4 
-0.9 


130-145 


1 


136.5 


0.1090 


+4.2 

-4.3 


±0.3 


+3.1 
-3.4 


+2.3 

-2.4 


+0.9 
-0.9 


+ 1.5 
-0.9 


50-145 


22 


M z 


0.1161 


+4.1 
-4.8 


±0.1 


+3.4 
-3.3 


+ 1.0 

-1.6 


+ 1.1 
-1.2 


+2.5 
-2.9 



contributions are estimated from the [i r j dependence of 
the calculations when varying the scales in the range 
0.5 < fi r .f/pT < 2. In the kinematic region under study, 
variations of fi r and /i / have positively correlated effects 
on the jet cross sections. A correlated variation of both 
scales is therefore a conservative estimate of the corre- 
sponding uncertainty. Since the [i r j uncertainties can 
not be treated as Gaussian, these are not included in the 
Hessian % 2 definition. Following Refs. [U, HH, the a s 
fits are repeated for different choices (fi r j — 0.5pr and 
Hrj = 2px) and the differences to the central result (ob- 
tained for [i r j = pt) are taken to be the corresponding 
uncertainties for a s (Mz)- Those are added in quadrature 
to the other uncertainties to obtain the total uncertainty. 

Data points from different \y\ regions with similar px 
are grouped to determine the results for a s (Mz) and 
&s(pt)- A combined fit to all 22 data points yields 
a s {M z ) = 0.1161+^8 with X 2 / N df = 17.2/21. The re- 
sults are shown in Fig.[T]as nine a s (pT) (top) and a s (Mz) 
values (bottom) in the range 50 < pr < 145 GeV with 
their total uncertainties which are largely correlated be- 
tween the points. Also included are results at lower px 
from inclusive jet cross sections in DIS from the HERA 
experiments HI HH and ZEUS [H] and the 3-loop RGE 
prediction for our combined a s {Mz) result. Our o; s (pt) 
results are consistent with the energy dependence pre- 
dicted by the RGE and extend the HERA results to- 
wards higher pt- The combined result is consistent with 
the result of a s (M z ) = 0.1189 ± 0.0032 from combined 
HERA jet data [3j] and with the world average value of 
a s (M z ) = 0.1184 ± 0.0007 ||. The contributions from 
individual uncertainty sources are listed in Table H The 
largest source is the experimental correlated uncertainty 
for which the dominant contributions are from the jet 
energy calibration, the pt resolution and the integrated 
luminosity. 

Varying the size of the uncertainties of the nonper- 
turbative corrections between a factor of 0.5 and 2 



changes the central value by _q odio an< ^ does not af- 
fect the uncertainty of the combined a s (Mz) result. Re- 
placing the MSTW2008 NNLO PDFs by the CTEQ6.6 
PDFs changes the central result by only ±0.5% which 
is much less than the PDF uncertainty. Excluding the 
2-loop contributions from threshold corrections and us- 
ing pure NLO pQCD (together with MSTW2008 NLO 
PDFs and the 2-loop RGE) gives a result of a s (M z ) = 
0.1202lQ QQgg. The small increase in the central value 
is a result of the missing 0(a A s ) contributions which are 
compensated by a corresponding increase in a s . The dif- 
ference to the central result is well within the scale un- 
certainty of the NLO result. The increased uncertainty 
is mainly caused by the increased /j, r j dependence, but 
also by the larger PDF uncertainty at NLO. 

In summary, we have determined the strong coupling 
constant from the inclusive jet cross section using the- 
ory prediction in NLO plus 2-loop threshold corrections. 
The o: s (pt) results support the energy dependence pre- 
dicted by the renormalization group equation. The com- 
bined result from 22 selected data points is a s (Mz) — 
0.11611q oo48- This is the most precise a s result obtained 
at a hadron collider. 
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